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Abstract 



We present explicit expressions for the correlation functions of interacting 
fermions in one dimension which are valid for arbitrary system sizes and 
temperatures. The result applies to a number of very different strongly corre- 
lated systems, including mesoscopic quantum wires, quantum Hall edges, spin 
chains and quasi-one-dimensional metals. It is for example possible to calcu- 
late Coulomb blockade oscillations from our expression and determine their 
dependence on interaction strength and temperature. Numerical simulations 
show excellent agreement with the analytical results. 
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Recently, there has been great interest in strongly correlated systems of all kinds, which 
is spurred by high temperature superconductivity, progress in mesoscopic physics, quantum 
Hall systems, and newly available materials which accurately display the characteristics of 
one-dimensional metals or spin chains. In the case of quasi-one-dimensional electron systems 
much progress can be made with the Luttinger liquid formalism |I| which is able to describe 
virtually any type of local interaction and can be solved by use of bosonization techniques 
Using these methods, it is well known how to calculate correlation functions in the 
infinite length or zero temperature limit and therefore determine the spectral properties 
Q of a number of systems, like mesoscopic quantum wires, quantum Hall bars |J, quasi 
one-dimensional metals and spin chains ||[7| even in the presence of boundaries HQ . We 
now present explicit expressions for the correlation functions which are valid at arbitrary 
temperature, system size, and distances (as long as the lattice spacing is relatively small). 
These expressions can be used to calculate Coulomb blockade oscillations in mesoscopic sys- 
tems as a function of the interaction strength and of temperature. Monte Carlo simulations 
of one of the simplest systems, namely interacting spinless fermions on a one-dimensional 
lattice, show excellent agreement with our result. 

As our model Hamiltonian we consider one-dimensional interacting fermions in the con- 
tinuum limit 

,/ 

H = dx 
Jo 

Here Jl/r = : ^l/r^l/r '■ are the chiral fermion currents of the left- and right-moving 
components of the fermion field *&(x), expanded about the Fermi points dbfep, = 
e- ikFX ip L {x) +e ikFX ip R {x). This expansion is valid as long as the lattice spacing a is small 
compared to the inverse temperature v/3 and the length scales we want to consider. In 
this limit the Hamiltonian (|l|) can describe almost any one-dimensional fermion system with 
short range interactions by choosing appropriate coupling constants gi and g 2 - The umklapp 
process e likpX ^l^r^liPr + h.c. can also be included, but it only contributes if the system 
is close to half-filling kp = it /2a and then it can often be absorbed by renormalizing the 
coupling constant g 2 0]. 

The fermion correlation functions of systems with the Hamiltonian of equation ([!]) are 
known to obey simple power-laws in the zero-temperature and infinite length limit, with 
exponents determined by the interaction 0. At non-zero temperatures the correlation func- 
tions are exponentially damped and are described by powers of the hyperbolic sine in the 
infinite length limit. For finite length I and periodic boundary conditions, on the other hand, 
the correlation functions are periodic and at zero temperature are described by powers of 
the sine instead. The correlation functions are essential to understand the spectral charac- 
teristics and other properties of experimental systems, and it is often necessary to consider 
both finite size and finite temperature. We will now derive explicit expressions for this case 
and show the crossover between the two limits. 

To establish notation we will review shortly the bosonization formalism H for the Hamil- 
tonian (HD, which is given by the relations 

Jl,r = ^~d x {ct> ± }), ip L>R oc e -*VSF(*±*) ( 2 ) 



d d 
i-T^L - v f WRi-j-^R + 9i{JlJl + JrJr) + 92JlJr 

\JjJU \AjJu 
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where <fi is the dual field to the boson <fi with the finite length mode expansion 



1 



Vt 



{x,t) = 4> + + + J2 



n =i 



e-¥«(«*+*) a ^ + e -^ n(vt ~^a^+ h.c. 



(3) 



Here II and Q are canonically conjugate to the zero modes 0o an d 4>o, respectively, and the 
time variable implicitly carries a small ultraviolet cutoff t — ia. 

The gx interaction can be absorbed by redefining the velocity v = Vp + gi/2ir, and the 
Hamiltonian density becomes 



n 



+ 



92_ 
47T 



- [(<9,0) 2 + (d x <P) 2 
The Hamiltonian is then solved by a simple rescaling of the boson, 



(4) 



(5) 



where K = 1 — gz/Anv to lowest order in the coupling constant (i.e. K = 1 for a non- 
interacting system and < 1 for repulsive interactions). 

The eigenvalues of the operators n and Q are quantized by using the relation for periodic 
boundary conditions t/>_r,l(0, t) = e ±zkFe ipR^L(^,t) together with equations (§) and @. After 
a careful consideration of the commutation relations we find 



Q = ^( n -n )/K, U = ^ttK 



m, 



(6) 



where n and m are either both even or both odd integers and no = — + 1 which can 
be defined modulo 2. The quantum number n represents the particle number and m is a 
measure of the current in the ring. If a magnetic field is applied, the ground state may have 



quantum number m^O and the system carries a persistent current [|10|,|IT 



As a simple example of our calculation at finite temperature and size, we now want to 
consider the chiral Green's function 



G(x,t) = (M^t)^ L (0,0) 



(7) 



but the generalization to more complicated correlation functions is straightforward as we 
will see later. Using equations (Q) and (|5|) and applying the Baker-Hausdorf formula as well 
as the cumulant theorem for bosonic modes, we obtain 



G(x,t) oc (exp -iy/Tr\K<f>(x,t)+}(x,t)/K] expi^ \K<j>(0, 0) + 0(0, 0)/K 



(8) 



-ij[x+- 



ot] B(x,t) exp\ir(K + K 1 ) 2 Gl(x, t) exp \tt(K — K l ) 2 Gn(x,t) 



Kvt+x/K 



where B(x, t) = ^exp —i\pn (O Kx+ J t / K _|_ n 
tors n and Q, and Gl,r = \4>L,R{x,t)<f>L,R{Q,fy — 0|,r(O,O)) are the chiral boson Green's 



is the contribution from the opera- 
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functions for the left- and right-moving modes 4>l,r = [(f) ± - (Q ± n)2f*]/2, respectively. 
Using the mode expansion (H) the chiral boson Green's function can be expressed as 



G L (x,t) 



4 7 r n e i nL 



e"**™^ - 1) (afc^) + (e i2 ™^ 



(9) 



Since we know the energy levels E n = n^f- of the free boson modes, it is useful to write 
the Bose-Einstein distribution as a sum over occupation numbers (a) n a r \ - 1 



YHk=i e ™ ^ kn . After exchanging the order of the summations over k and n we can use 
the identity E^Li v = - hi(l - x) to obtain 



G L {x,t) = -—In 

47T 



]_ _ e -27ri(x+^)/£ oo -y _ 



-27TJ 



(x+vt-iv/3k)/e y _ e 2m{x+vt+ivl3k)/l 



I — e -2nv(3k/e 



I — e -2iTvf3k/l 



(10) 



Using the elliptic theta function of the first kind 9%(z, q) [|T^], this result can be written more 
compactly as 



G L {x,t) 



.X + Vt 1 ^i(7T , . 

', In £ — 



X+Vt g— ■KV 



(11) 



and likewise for Gn(x,t) = Gl(— x, t). 

The cumulant theorem does not apply to the contribution B(x,t) from the operators II 
and Q, but we can use the energy spectrum E = %i(Q 2 + II 2 ) and sum over all eigenvalues 
in equation (H). The result can again be expressed in terms of the elliptic theta functions 





B(x,t) 



JnoXQ e 2 (x Q +i^, e-^ K2 )6 2 (xu, e- 2 ^ 2 ) + e 3 (x Q + i^, e-^ K2 )6 3 (xu, e- 2 ^ 2 ) 



e 2 {ilfi, e- 2 -r/ K2 )6 2 (0, e- 2 ^ 2 ) + 3 (i7ffe, e - 2 ^ 2 )# 3 (0, e" 2 ^ 2 ) 



(12) 



where 7 = is a measure of the energy gap compared to the temperature scale, and 
we have defined xq = j(x + vt/K 2 ) and in = j(x + vtK 2 ). The complete chiral Green's 
function therefore reduces to a compact expression 



G(x, t) oc B(x, t) 



0i (tt 



vt—x p — nv 



l-(K--k) 2 /4 



(13) 



We can immediately verify that this expression is anti-periodic under translation t — > t + i/3 
and periodic under translation x — > x + £ up to a phase of e* fcF£ as it should be. 

If the temperature is much smaller than the energy gap of the system, i.e. 7 = nf- -> 00, 
we can use the limits of the theta function [12 as q = e -7 — ► to find 



G(x,t) 



3 i(no mod2)a;Q 



t . . x + vt, 
-sin vr — — 

7T £ 



n-(^+^) 2 /4 



vt — x s 



— sml 7T- 

7T 



-(i^-i) 2 /4 



(14) 
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which is the expected finite length result [0J. For special values of no the overall phase may 
have a different x-dependence, which will not be discussed here in detail. 

Likewise, we can explore the limit of a large system size compared to the tempera- 



ture scale, 
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i.e. 

2,f^e- z2 ^e- w2 ^smh^ and 6 2 

V 7 7 

correct finite temperature result 



0. We use the Poisson summation formula to find 9\ 

t) pa 0of* ~ pKp-^h 



z,e ') w y 3 ^,e 
fl3|l is reproduced 



-e 

7 



as e 



z, e 'J pa 
so that the 



G(z,t) 



u/3 x + ft , 
— smh(7r — 

7T ^ Vf3 ' 



-(K+i) 2 /4 



v/3 vt — x 
— smh vr — 

7T Up 



(15) 



Taking the limit I — > oo in equation fll4]) or u/3 — > oo in equation fll5|) recovers the well 
known universal power-laws. 

One interesting aspect of equations 0121) and (jl3f) is the dependence on the variable 
no = ~ E -+1. The Green's function (|T3|) is invariant under the increase of no by 2, i.e. periodic 
in kp with period Hence, as we change the Fermi energy (e.g. by applying a gate 
voltage), the spectral properties of the system are changed periodically. These are Coulomb 
blockade oscillations 



15,16 



which come from the fact that the system has a regular energy- 
level spacing. In the case of Luttinger liquids this level spacing has been derived from a 
complete solution of a quantum system which takes all interactions into account and cannot 
be obtained from a geometrical analysis of the capacitance as in the usual quantum-dot 
systems. As possible experimental setups to test this periodicity we can imagine electron 
tunneling through a quasi one-dimensional wire or ring of mesoscopic size. Such wires can 
be fabricated by etching, gating [T7| or by cleaved edge overgrowth flIB| . A small gated 
quantum hall bar on the other hand should exhibit chiral Luttinger liquid behavior which 
can be examined with the same formalism, except that in that case the rescaling parameter 
K is fixed by the filling fraction in the bulk of the quantum Hall bar ||. The experimental 
situation is developing quickly and most recently carbon nanonubes have been produced 



which were seen to exhibit the charactersitcs of one dimensional wires 

As an example we consider the Green's function at the lowest Matsubara frequency 
iuo = in/ (3 for spin-less Fermions on a ring. This quantity only gives an indirect indication 
of experimental tunneling resonances, but it is a simple illustration of the properties of 
equation ( Jl3| ) as a function of temperature and interaction. The absolute value of the 
Matsubara Green's function has been plotted in figure (P as a function of the Fermi wave- 
vector in arbitrary units. As the Fermi-level is changed, the system shows periodic resonances 
which become more pronounced as the temperature is lowered. These peaks occur at special 
values of no = 1+ 2 X , where the even and odd sectors of n and m in equation (^) give 

the same degenerate ground state. Hence, the spacing of the resonances can give a direct 
experimental measure of the interaction strength K [|lTJ. As shown in the inset of figure ([!]) 
a non-interacting system (K = 1) has only one central peak in the range < n < 2, while 
an interacting system has two split resonances. The physical interpretation of the resonances 
is that tunneling of an extra electron into the system does not require any additional energy 
at the degenerate values of n so that such processes are enhanced. Recently, much progress 
has also been made with Bethe ansatz methods |T^], which showed evidence of Coulomb- 
blockade-like oscillations in a strongly interacting system, but in contrast to statements 
made in that paper, field theoretical methods also predict the desired Coulomb blockade 
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oscillations fll5l , |16| . Using our result (|13D it is now possible to determine the complete spectral 



properties, including the dependence of the Coulomb blockade oscillations on interactions 
and on temperature for any of the systems mentioned above. 

So far we have taken the fermion field to be spinless, but the results apply equally well to 
electrons with spin because the well known spin-charge separation allows the same formalism 
to be used for the spin and the charge excitations separately [Q. It is also well understood 
how to incorporate open boundary condition within the same formalism P,PU[ • 

To test our result, we will now consider one of the simplest models of spinless interacting 
fermions described by the Hamiltonian 



H = - 



tJ2 [p\x)i)(x + 1) + h.c] +uJ2 



n[x) 



n{x + l) - - 



(16) 



where n = ip^ip, and the lattice constant has been set to unity. In this case the coupling 
constants are given by g\ = g-ijk = U. We have chosen the system to be at half-filling 
kp = 7r/2 where the model is equivalent to the xxz spin chain. In this case the umklapp 
process is allowed, but irrelevant, and the rescaling parameter K and the velocity v can 
actually be determined exactly by comparison with Bethe ansatz results ||21|| , i.e. to all 
orders in the coupling constants. For our numerical simulations we have taken t = U 
which corresponds to K 2 = | and v = i^rp. In particular, we consider the equal time 
density-density correlation function (n(x)n(O)) . At half filling the density- density correlation 
function acquires an alternating part which dominates the uniform contribution and can be 
expressed in terms of chiral fermions 



(n(x)n(O)) 



alt 



oc e 



2ik F x / / 1 



(17) 



Following the analogous steps from equation (§D to equation (|13|) we find 



(n(x)n(O)) 



alt 



OC 



0i(¥,e- 



#i(T> e " 7 ) 



-3/2 



(Mjf, e- 8 ^ 3 )fl 3 (0, e- 3 ^ 2 ) + 3 (gf , e- 8 ^ 3 )9 2 (0, e^/ 2 ) 
2 (O, e- 8 ^ 3 )^ 3 (0, e-^/ 2 ) + 3 (O, e- 8 V 3 )e 2 (0, e" 3 ^ 2 ) ' 

(18) 



Here we took n to be odd, i.e. £ divisible by four, but other cases have similar expressions. 

For comparison we performed numerical Monte Carlo simulations using the Hamiltonian 
(|IT)|) and extracted the alternating part of the density- density correlation function to de- 
termine the unknown proportionality constant in equation (|18"D. This prefactor should be 
independent of x up to corrections for smaller distances, which come from irrelevant terms 
in the Hamiltonian that can be neglected in the long-distance limit x 3> a. As shown in 
figure @ the prefactor is indeed constant within the error-bars of the Monte Carlo simu- 
lations for all system sizes I and temperatures. This serves as a very sensitive test of our 
theoretical predictions since the expression fll8|) varies over several orders of magnitude over 
the range shown, so that even slight changes from the predicted expression flI8|) would have 
resulted in huge deviations from a constant for larger values of x. Using the previously 
known approximate expressions (fH| ) or ( Jl5| ) is not sufficient in this case. 

We therefore conclude that our field theoretical calculations make very accurate estimates 
for the correlation functions of almost any one- dimensional interacting system at arbitrary 
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temperature and system size (as long as the cutoff a is small compared to the corresponding 
scales). The results can be used to determine the complete spectral behavior of such systems, 
including Coulomb blockade oscillations as a function of both interaction and temperature. 
Monte Carlo simulations serve as accurate "experiments" , which confirm the results. 
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FIGURES 
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FIG. 1. The absolute value of the Green's function at the Matsubara frequency iujo = m/f3 
as a function of the Fermi wave- vector /cp for different temperatures and K 2 = 0.7. The periodic 
resonances in the spectral properties are the so-called Coulomb blockade oscillations. The inset 
shows the curves for different interaction parameters K 2 = 0.5, 0.7, 1 at ^fi = 30. 
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FIG. 2. The prefactor of equation fll8|), which is determined by dividing the alternating part 
of the density-density correlations from numerical simulations by the expressions of equation (18) 
for each correlation length x separately. As predicted, the prefactor is asymptotically constant in 
all cases. 
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